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Abstract 

We present and constrain a cosmological model where the only com- 
ponent is a pressureless fluid with bulk viscosity as an explanation for the 
present accelerated expansion of the universe. We study the particular 
model of a bulk viscosity coefficient proportional to the Hubble parameter: 
Cm = C^^) where C ~ constant. The possible values of are constrained 
using the SNe la Gold 2006 sample, the CMB shift parameter R from the 
three-year WMAP observations, the Baryon Acoustic Oscillation (BAO) 
peak A from the Sloan Digital Sky Survey (SDSS) and the Second Law 
of Thermodynamics (SLT). It was found that this model is in agreement 
with the SLT using only the SNe la test. However when the model is con- 
strained using the three cosmological tests together (SNe-l-CMB-l-BAO) 
the results are: 1.- The model violates the SLT, 2.- It predicts a value of 
Ho ~ 53 km • sec~^ ■ Mpc~^ for the Hubble constant, and 3.- We obtain a 
bad fit to data with a Xmin ~ 532 (Xd.o.f. ~ 2.92). These results indicate 
that this model is viable just if the bulk viscosity is triggered in recent 
times. 

We present a flat cosmological model which component is a pressureless fluid 
made of baryon and dark matter components with a bulk viscosity coefficient 
^ni proportional to the Hubble parameter, i.e.. Cm = CH, where C = constant 
and H is the Hubble parameter, to explain the present accelerated expansion 
of the universe. The subscript "m" stands for matter with bulk viscosity. A 
bulk viscosity coefficient can produce a positive term in the second Friedmann 
equation that induces an acceleration [1] [5] (i.e., a > , with a denoting the 
scale factor and the dots derivatives with respect to the cosmic time). Similar 
models or analysis have been proposed also in [51 [31 SI O [S] . 

A bulk viscous matter-dominated universe The energy-momentum ten- 
sor of a fluid composed by only matter (hereinafter we call matter to baryon 



and dark matter components together) with bulk viscosity Cm is defined by 
Tf_,tj = PmUf_,u^ + [g^u + u^Ui,)P^, where P„* = - 3Cm (a/a) (see e.g. [1]). 
Here Pm is the pressure of the matter fluid and the function a the scale fac- 
tor. The four-velocity vector is that of a comoving observer that measures 
the pressure and the density of the matter fluid. Assuming pressureless 
matter Pm = imphes P^ = — 3Cm {a /a). 

On the other hand, the conservation of energy has the form a{dpm/da) = 
3 {SCaiH — pm), where H = d/a is the Hubble parameter. Using the ansatz 
Cm — CH for the parametrization of the bulk viscositj0 and the first Fried- 
mann equation — (87rG/3)p„i (assuming a spatially fiat geometry for the 
universe), the conservation equation becomes a{dpm/da) = Spm [(SttG/S) C ~ !]• 
The solution in terms of the redshiflU 'z' is Pmiz) = p'^{l + z)^^^ , where we have 
defined the dimensionless bulk viscous coefficient ( = SttGC and is the value 
of the density of matter evaluated today. We substitute this solution into the 
first Friedmann equation so it becomes = (SttG/S) p°-j(l -I- z)^~''. Dividing 
the last expression by the critical density today = S-ffg / SttG we obtain 

H{z)^Ho{l + z)^('-^^ (1) 

where Hq is the Hubble constant. In this model the matter density parameter 
is rij^ = Pm/Pcrit = ^ becausc the matter is the only component of the universe. 



Cosmological tests We constrain the possible values of the bulk viscosity 
C using the following cosmological tests: the Gold 2006 SNe la data sample 
[7] composed by 182 SNe la, the Cosmic Microwave Background (CMB) shift 
parameter 'P' from the three- year WMAP observations [5], the Baryon Acoustic 
Osciflation (BAO) peak 'v4' from the Sloan Digital Sky Survey (SDSS) ^9^ and 
the Second Law of Thermodynamics (SLT) [IJ. 

For the SNe la test it is defined the observational luminosity distance p^fTT] 
in a fiat cosmology as c?l = c{l + z)Hq^ E{z')~^ dz' , where E{z) = H[z)/Hq 
and c the speed of light. The theoretical distance moduli for the z-th supernovae 
with redshift Zi is p{zi) = 51og]^Q[c?i(zi)/Mpc] -I- 25 . The statistical function 

r ~ n 2 

XsNe becomes XsNe(C-f^o) = Efc=i p(^fe, C, -f^o) - Pfe /o-^, where pk is the 
observed distance moduli for the fc-th supernovae and cr^ is the variance of the 
measurement. 

The CMB shift parameter R is defined as P = ^^0^// P(z')~^ dz', 
where ^cmb — 1089 is the redshift of recombination [51[T^. The observed value 
of the shift parameter R is reported to be Pobs = 1-70 ± 0.03 [S]. 

From the SDSS data, we use the baryon acoustic oscillation (BAO) peak A 

defined as A = ^/9^^E{zi)-^/^ [z:['^ E{z')-^ dz'f^ , where zi = 0.35 0. 
The observed value of SDSS-BAO peak is ^obs = 0.469 ± 0.017 [9J. 



^The explicit form of the bulk viscosity is assumed a prior. See for instance [2l [3l l5l |6) . 
^The superscript "0" (zero) represents a quantity evaluated today and the relation between 
the redshift z and the scale factor ais(l-|-2) = l/a. 
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Figure 1: Confidence intervals for the parameters (CHq) for a spatially flat 
bulk viscous matter-dominated universe. We show the confidence intervals of 
68.3%(1ct), 95.4%(3cr) and 99.73%(5cr). The constraints are derived from the 
Gold 2006 SNe la data sample alone (left), and the joint Gold 2006 + CMB + 
BAO cosmological tests (right). Hq is in units of km • sec~^ • Mpc~^ and ( is 
dimensionless {( = 8TrG(). 



We construct the total function as xLtai = Xsnc + Xcmb + XbaO' where 
XsNe was defined above, Xcmb = [(^ - Rohs)/crR? and xIao = " 

The law of generation of local entropy in the space-time is found to be 
TWys" = C,m^uu'' ~ 5H(^ where T is the temperature, Vi^s" is the rate at 
which entropy is being generated in a unit volume. The second law of thermo- 
dynamics can be written as TV^s'^ > which implies that 3i?Cm > 0. Since H 
is positive in an expanding universe then = has to be positive in order to 
preserve the validity of the second law of thermodynamics. Thus, in the present 
model a necessary condition is C > 0. 



Constraining the bulk viscous cosmological model We compute the best 
estimated values and "the goodness-of-fit" of C and Hq to the data through x^- 
minimization, using SNe la test alone, and also considering the joint SNe la, 
CMB and BAO tests together. Then we compute the confidence intervals for 
{(,Hq) to constrain their possible values. We obtain as best estimates using 
only the SNe la test: C 1-43 ± 0.12 and Hq = 61.758 ± 0.86, with a x^n = 
165.135 (Xd o f = 0.917), where Hq is in units of km • sec~^ • Mpc^^. On 
the other hand, using the SNe la, CMB and BAO tests together we obtain 
C = -0.1451±0.02 and Ha = 53.425±0.4, with a xLn = 532.056 (xlo.f. = 2.92). 
It can be seen that when we consider only the SNe la test we obtain C > with 
a 99.7% confidence level and a reasonable value for the Hubble constant in the 
interval 58.9 < Hq < 64.8 as well as for Xd o f ~ 0.917. However, when the same 
analysis is performed using the three cosmological tests together (SNe la, CMB 
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Figure 2: Confidence intervals for the parameters (C, Hq) for a spatially flat 
bulk viscous matter-dominated universe. We show the confidence intervals of 
68.3%(lcr), 95.4%(3cr) and 99.73%(5cr). The constraints are derived from the 
Gold 2006 SNe la data sample alone and the joint Gold 2006 + CMB + BAG 
cosmological tests. Hq is in units of km • sec^"'^ • Mpc^^ and C, is dimensionless 

(C^s^GC). 

and BAG) we obtain negative values of C with at least 99.7 % confidence level 
that disagree with the second law of thermodynamics, a not so good estimation 
for Ho of 52.1 < Hq < 54.8, and a bad Xd o f = 2-92. These results are 
illustrated in figures [T] and [2] 



Therefore, we conclude of this model that: 

• In order to have a viable model for explaining the observed accelerated 
expansion of the universe, the bulk viscosity of the fluid with the ansatz 
considered in the present work should be triggered in recent times (z < 2). 

• Since the GMB and BAG tests have information for very large redshifts 
{z ^ 2) and due to the fact of that this model has a bad fit and violate 
the second law of thermodynamics when we test it using the joint SNe 
la, CMB and BAG tests then we conclude that this is not viable if it is 
extrapolated for early times of the universe as suggested in 0. 
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